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Abstract 



Correlation functions of the XXZ spin chain in the critical regime is studied at zero- 
temperature. They are exactly represented in the Fredholm determinant form and are related 
£N) , with an operator-valued Riemann-Hilbert problem. Analyzing this problem we prove that a 

C^l ' two-point correlation function consisting of sufficiently separated spin operators is expressed 

by power- functions of the distance between those operators. 

o 

^ ! 1 Introduction 

B 

1 ■ We consider a two-point correlation function of the XXZ spin chain in an external magnetic 

field at zero-temperature. The Hamiltonian is defined by 

8 : h xxz = Y,(«+i + «+i + - h<)- (1.1) 

The anisotropy parameter A is parameterized as A = cos 2?? < r/ < ir), which implies the 
critical regime. The Pauli matrices a^, cr^ act on the n-th site and h indicates an external 
magnetic field. The XXZ spin chain is one of quantum solvable models and is analyzed by means 
of Quantum Inverse Scattering Method Based on this method the generating functional of 
correlation functions is represented in the Fredholm determinant form j2|. 

The determinant representation of the generating functional has been investigated in some 
limiting cases. In the limit of strong magnetic field h — > h c = 4 cos 2 r\ this is connected with the r- 
function of the Painleve V equation || . On the basis of this equation we can evaluate correlation 
functions of the XXZ spin chain in strong magnetic field 0]. The generating functional depends 



on an extra parameter a (see (2.2)) and in the limit a — > — 00 itself indicates a special correlation 
function called the ferromagnetic-string- formation-probability (FSFP). The FSFP describes the 
probability of finding a ferromagnetic string of adjacent parallel spins and the corresponding 
Fredholm determinant is known to be related with an operator- valued Riemann-Hilbert problem 
d]. Thus the evaluation of FSFP is reduced to the solution of this problem. We remark that 
either the analysis of the Painleve V equation or Riemann-Hilbert problem is a kind of classical 
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inverse scattering problems. For special cases above, the calculation of correlation functions of 
the XXZ spin chain can be reduced to a classical inverse problem, despite the model is quantum 
mechanical. 

Recently, the operator-valued Riemann-Hilbert problem associated with the generating func- 
tional has been proposed ||. Surprisingly, the calculation of any correlation function of the 
XXZ spin chain is interpreted as a classical inverse scattering problem. In particular we are 
in a position to obtain the closed forms of two-point correlation functions: there was no valid 
method to evaluate two-point correlation functions of the XXZ spin chain from the stand- 
point of quantum solvable models, because the integral representations of two-point correlation 
functions obtained in |?J and Hj are too complicated to analyze them. In this paper, on the 
basis of our Riemann-Hilbert problem, we compute the large-m asymptotic form of a two-point 
correlation function (cr^crf) and prove the following theorem: 

Theorem. For sufficiently large m a two-point correlation function (a^af) is expressed by 
power-functions of m. 

This corresponds to a rigorous proof of the well-known conjecture that two-point functions of 
the XXZ spin chain in the critical regime decay in power of the distance between spin operators 
at zero-temperature, which was derived by means of Conformal Field Theory, the bosonization 
method, the numerical analysis and so on. In order to prove the theorem, in section 0, we 
define the Fredholm determinant representation of the generating functional and formulate the 
associated operator- valued Riemann-Hilbert problem. In section |3] we analyze this problem and 
determine the asymptotic form of generating functional. As a result the theorem is proved. The 
last section ^| is devoted to concluding remarks. 



2 Operator-valued Riemann-Hilbert problem 

In this paper we pay attention to a two-point correlation function (a^af), which is computed 
through the generating functional. In this section we define the Fredholm determinant repre- 
sentation of the generating functional and formulate the associated operator-valued Riemann- 
Hilbert problem. 



2.1 Fredholm determinant 

A two-point correlation function (a^af) is given by the generating functional Q^ m \a) via [Q 



(a z m al) = 2A r 



d 2 Q( m \a) 



da 2 



9Q (1) ( 



a) 



a=0 



da 



+ 1. 



(2.1) 



Here A m is the lattice Laplacian defined by A m f(m) = f(m) — 2f(m — 1) + f(m — 2). The 
generating functional Q^ m \a) is represented by the Fredholm determinant: 

Q( m )(a) = Q (vacjDet(l - ^( m ))|vac) (2.2) 

where Qq is a normalization factor. The Fredholm determinant is expanded as 



00 1 / 1 \ n r n 

Det(l-£VM) = X: -J --) /n^det, 

n=0 ' v ' k=l 



(2.3) 



2 



where det n V^ is the normal determinant of an n x n matrix Zj , Zj 



(i,3 
-AO 



1, 



,n). The 

integration contour T runs on a part of the unit circle anti-clockwise: T B e iy (2tt — ip<9<ip). 
The boundary value tp depends on the anisotropy parameter 77 and the external magnetic field 
h, lies in it < tp < 2rj and approaches to tt in the strong magnetic field limit h — > h c = 4 cos 2 rj, 
to 2r] in the limit h -> [|, ||. The kernel is given by the scalar product of four-component 
vectors as follows Hi: 



V {m \z, 
where vectors are defined by 



w 



ds 



^4 m) (z| S )^ m) ( W | S ) 



w 



(2.4) 



a 



(m) 



(z\s) 



1/2 



z m / 2 exp(ig- 1 ^s + a- 



z 



z 

-m/2 



Ps(*)) 



exp^" 1 ^*-^)) 



V z- m / 2 exp(-iq^-s + a + Mz)-^3(z)) J 



z\s 



i(fl 



1/2 



V 



^ m/2 exp(-i^ s + ^ 4 (z)) 

exp(-i^ S + ^ 2 (z)) 

r m / 2 exp^- 1 ^^- 



(2.5) 
(2.6) 



and q = e 2 " 7 . The operators </?j(-z) (i = 1, • • • ,4) are bosonic quantum fields called the dual fields. 
They are decomposed into coordinate and momentum fields as <fj(z) = qj(z) + Pj(z) and act 
on the vacuum states defined by (v&c\qj(z) = pj(z)\v&c) = (j = 1, ... ,4). The commutation 
relation between coordinate and momentum fields is expressed in the following matrix form: 



[qj(z),p k (w)] 



( 1 




V 1 



1 \ 



1 
1 1 



log h(z, w) + 



with 



h(z, w) 



11/ 



1 



( 1 



1 



1 \ 



1 1 

1 

1 1 




1 

1 / 



log h(w, z) 



(2.7) 



q — q 1 V~V z — q qw — 1 \ qz — I w — q 

By definition the dual fields are commutative: [<Pj{z), <pk(w)] = (j, k = 1,...,4). They are 
thus interpreted as entire functions unless they act on the vacuum states to produce the vacuum 
expectation values. The following lemma is useful. 

Lemma 1. A dual field 924(2) is equal to ^3(2). 



qz — 1 w 



q qw — 1 



2i 



Proof. If a set of bosonic fields with four-species has the same commutation relation as ( |2.7| ), 
the kernel whose dual fields are replaced with those fields is equivalent to the original 

kernel. We put c-numbers a and £ ±m / 2 into the dual fields. Under a simultaneous replacement 
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<Pl P2 and y>3 <-> y?4 the commutation relation (p.7| ) is invariant and the kernel W" 1 ) (z, io) is 
then changed into exp(<^3(w;) — (p±(z))V( m > (w, z). It thus follows that 

V^(z,z)=ex V ( V3 (z)-M^))V (m) (z,z). (2.9) 

Note that the kernel V^ m \z, z) contains the dual fields <pi(z), . . . , (f4(z). In order to produce the 
same vacuum expectation values in both sides, a dual field 924(2) must be equal to fz{z). □ 

Hearafter we replace a dual field ^4(2) with <f3(z) on the basis of this lemma. By virtue of 
this lemma the recursion relation (|3.21|) contains no dual fields and makes clear that two-point 
correlation functions decay in power. 

2.2 Operator- valued Riemann— Hilbert problem 

The Fredholm determinant Det(l — ^V^ 1 ) is related with matrix integral operators ||. They 
are defined by matrices whose multiplication requires extra integral calculus: those products are 
formulated as 

poo 

{EF)ij{z\8,t) = drJ2E ik (z\s,r)F kj (z\r,t). (2.10) 
Jo k 

From now on we omit the dependency on extra variables s, t of matrix integral operators unless 
necessary and denote the delta- function S(s — t) by I. Consider a 4 x 4 matrix integral operator 
X^ m \ z ) that satisfies the following conditions (a)— (c) Q: 

(a) Each element of x^ m \ z \ s ^) ls an analytic function of z G C\r for any m G Z>o, s,t G C. 

(b) x (m) (°o) = L 

(c) Let x+ \ z ) (X—( z )) he the boundary value of X (m) {z) on the left (the right) side of a 
positive direction of the oriented contour T. It then follows that 

X (m \z)= X ^\z)L^(z). (z G T) (2.11) 

The 4x4 matrix integral operator L^ 71 ' (z) is called the conjugation matrix and is expressed by 

/ I-kP(z) Ke-^Q(z) -Kz m e^-^P(z) Kz m e-^Q(z) \ 

-eWR(z) I + P T {z) -z m e^R{z) z m e-^ 1+ ^P T (z] 

-z- m e-^ + ^P{z) z- m e~ V2 Q(z) I — P(z) e - Vl ^ 2+V3 Q(z) 

\ -KZ- m e^R(z) Kz- m e^-^P T (z) -Ke^ + ^-^ 3 R{z) I + kP t {z) j 

(2.12) 

where k = e a and the z-dependency of the dual fields is omitted. The integral operators P(z), 
Q(z), R(z) and P T (z) are given below: 

P(z) = i(q-q~ 1 ) — exp ( iq^ " 



L^ m \z) 



Q(z) = i{q-q )— r exp(ig 

V z Q 

- 1 ■/ - ■ ( qz — 1 
R(z) = i{q-q ) exp [ -iq 

V z-q 

,->/', , •/ -i\ q* — x f . qz - 1 
P (z) = i[q — q ) exp I —iq 



qz 


- 1 


z 


- Q 


qz 


- 1 


z 


- Q 


qz 


- 1 


z 


- Q 


qz 


- 1 



z-q 
x qz-\ 



-s — iq 1 


(2.13) 


z-q J 


: (* + *)) 


(2.14) 


(* + *)) 


(2.15) 


s + ig" 1 ^— U ) . 


(2.16) 


z-g J 
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They satisfy the following relations: 



P{z)P(z) = P(z) P(z)Q(z) = Q(z) (2.17) 

Q(z)R(z) = P(z) Q(z)P T (z) = Q(z) (2.18) 

R{z)P(z) = R{z) R(z)Q(z) = P T {z) (2.19) 

P T {z)R{z) = R{z) P T {z)P T {z) = P T {z). (2.20) 

Here these products implicitly contain the integrals for extra variables s, t like for ( 2.10[ ). 

To find the 4x4 matrix integral operator that satisfies (a)-(c) is referred as the operator- 
valued Riemann-Hilbert problem. By using its solution x ( z ) the Fredholm determinant 
Det(l - ^V {m) ) is computed via the following lemma |H|: 

Lemma 2. The Fredholm determinant Det(l — -^V^) obeys 

DetCl LT/( m +!h , 

De ; (1 -V(.»,) ) =*W"'(0)]-S) P.21) 

where 72 = 3(1 — a z ® 1). 

Proof. By Lemma 5 and (4.7) of || we have 

Det(l-iy(-+ 1 ))___ l ^ r .»rmi-i. 



Det(l - i7H) 



expTrlog( 7 i + 72 [x lmj (0)]- i 72 ) (2.22) 



where 71 = |(1 + a z <8> 1). Due to the orthogonality 7172 = and 7i = 71, 7 2 = 72 it follows 
that 

expTrlog(7 1 +7 2 [ X {m) (0)]- 1 7 2 ) = expTr(log 7l + log( 72 [x {m) (0)]- 1 7 2 )) 

= expTrlog( 72 [x (m) (0)]- 1 7 2 ). (2.23) 

Using the relation logDet = Trlog we obtain the lemma. □ 

On the basis of this lemma we derive the large-m asymptotic representation of a two-point cor- 
relation function (a^af) from the solution x i z ) °f operator-valued Riemann-Hilbert problem 
(a)-(c). 



3 Asymptotic form of two-point correlation function 



In this section we compute the large-m asymptotic behavior of the solution ( z ) °f operator- 
valued Riemann-Hilbert problem (a)— (c). By virtue of the relations (2.17)-( 2.20[ ) the conjuga- 
tion matrix L^ m \z) can be decomposed into the product of triangular matrix integral operators 
as follows: 

L^ m \z) = M (m \z)N (m ^(z) 



where 



(I 





Vo 



I 





-z m e^-^P{z) 
-\z m ^R{z) 

I 





^z m e-^Q(z) \ 



-<Pi-<pa+tf>3Q( z ) 
I 



(3.1) 



(3.2) 
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and 

/ I-i^Piz) \ 

1 ,__ V - 1 „T / V 



jvM(j5 



^ew>i2(jc) I+\P T {z) 

\ -Kz~ m e^i?(z) Kz- m e^-^ 3 P T (z) -Ke^+^-^i^z) I + kP t (z) J 

(3.3) 

Let Ti (r2) be a contour that runs from e to e 1 ^ inside (outside) the unit circle. Consider 
a matrix integral operator x ( z ) defined below, like for the case of the lattice sine-Gordon 
model §§. 

(i) x {m) {z) = X (m) (z) outside the region enclosed by the contours T\ and T2- 

(ii) x^ m \ z ) = X^ m \ z )M^ m \z) m the region enclosed by T and Fi. 

(iii) x^ m \ z ) = X^C^X-^ 2 )) -1 m the region enclosed by V and T2- 
By definition analytic for z £ C\(Fi U T2) and obeys 

xL m) (z) = xt\ z )M {m H z ) (2€rx) (3.4) 
= ^M (z)A rM (z) . (zer2) ( 3. 5) 

We are interested in the asymptotic behavior of the solution x^ m \ z ) f°r sufficiently large m. 
Notice that in the limit m — ► 00 the conjugation matrices M( m )(z) and 7V( m )(z) become block- 
diagonal in the vicinities of contours T\ and T2, respectively. Take this limit and consider the 
original conjugation matrix for x^ m \ z ) again. Then we have 

* M <=> - ( T 4) ) • (3 - 6 > 

The 2x2 matrix integral operators <&i(z), $2(2) satisfy 

*-(*) = ®t{z)G a {z) (zeT, a = 1,2) (3.7) 
and the corresponding conjugation matrices are expressed by 



Gl(z) " ( ^R(z) I + \P T {z^ ' (3 ' 8) 



r ^ ( J - P ^ e-^-^g (z) \ 

~ V -Ke^+w-V3i?(z) / + K P r (z) )■ 

The calculation of a two-point correlation function (cr^crf) requires the Fredholm determinant 
of the solution x ( z )- Taking the Fredholm determinants in both sides of ( |3.7| ) we obtain the 
scalar Riemann-Hilbert problems: 

Det$~(z) = Det$+(z)DetG a (z). (z G V, a= 1,2) (3.10) 

In the same way as Q the Fredholm determinants of Gi(z), G2(z) are computed as 

DetGi(z) = k" 1 DetG 2 (z) = k. (3.11) 
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According to the condition at the infinity: Det<l?i(oo) = Det3>2(oo) = 1, the solutions of scalar 
problems ( 3.10| ) are determined as follows: 

They are analytic for z 6 C\r. 

We evaluate the m-dependency of the solution x^ m \ z )- Based on the decomposition ( |3.1|) 
the conjugation matrix for sufficiently large m is expressed as follows: 

L [Z) - { z- m G 21 (z) G 2 {z) ) [6A6) 

where Gi 2 (z), G 2 \{z) are 2x2 matrix integral operators independent of m. Let the solution 
X {%) denote by 



X^'{z) X^'(z) 



Here X^\z), X^\z) decay sufficiently rapidly as m — > oo. By the relation (2.11) of the 
original Riemann-Hilbert problem, the 2x2 matrix integral operators x[ m \z), Xj™\z) obey 

[xl m \z)U = [X^\z)] + G 1 {z) + z^X^ {z)] + G 21 {z) . . , , 

[Xt\z)\- = [Xi m \z)] + G 2 (z) + Z ™[xi?\z)] + G 12 (z). { ] { - ' 

In terms of the relation the conjugation matrix G a (z) is given by [<3?+ (z)] _1 3>~(.z) (a = 1, 2). 
Applying them we have 

[xj m >(20]_[*r(*)] -1 = [x^( 2 )] + [$+(^)]- i +^[x 1 ^( 2 )] + G 21 (^)[$r(^)]- 1 (3.16) 

[^(z)]-!*^*)] -1 = [^^)] + [$+(^)]- 1 + ^[X 2 ^^)] + G' 12 ( 2 )[$ 2 -(^)]- 1 . (3.17) 

Here the variable z lies on the contour T. We remark that the inverse of & a ( z ) exists because 
Det^ a (z) is non-zero for z e T (a = 1,2) (see (gig)). By means of the Plemelj formulae (see 
||Io|| ) and the condition x[ m \oo) = (oo) = /, the 2x2 matrix integral operators x[ m \ 
X^ (z) are expressed below: 

X[ m \z) = (/-^^^^[^Hl+G^H^rW]- 1 )^^) (3.18) 
Af^(z) = (/-^^^[^Hl+^H^W]- 1 )^^). (3.19) 

A two-point correlation function (cr^af) is computed by using only [x (m) (0)] _1 - In the case 
z = the right sides yield some terms whose order is less than 1/m, because the contour T is 
a part of the unit circle and x[^\z), X^\z) decay sufficiently rapidly asm-> oo. Hence we 
obtain the large-m asymptotic behavior of the solution x (0) as follows: 

(ra) /^(OKHOfm- 1 )) \ 

X (Uj "V $ 2 (0)(/ + O(m- 1 )) J ■ i m>>i ) ^.2Uj 
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We are in a position to compute a two-point correlation function (cr^af). Applying the 
expressions fl3.12|) and ( gjjj) to Lemma | we have 



^ 2 \ . J =exp [Z-a) (l + Oim- 1 )). (m » 1) (3.21) 

Det(l-^F( m )) / 

Since the leading term exp(tpa/ir) does not contain any dual field, the vacuum expectation value 
of the left hand side can be written as 



Q( m+1 >(q) 
Q(™)(a) 



exp (^-a\ (l + 0(m" 1 )). (m > 1) (3.22) 



From this recursion relation we derive the large-m asymptotic form of the generating functional 
as follows: 

Q( m )( Q ) = f(m, a) exp f ^«™J • (m » 1) (3.23) 

Here f(m,a) satisfies /(m + l,a)//(m, a) = 1 + 0(m _1 ). Substituting this asymptotic form 
into the relation (|2.l[) we arrive at the following relation for two-point correlation functions: 

/ _2 _Z\ 

l + Otm- 1 ). (m>l) (3.24) 



Thus the asymptotic form of a two-point correlation function of the XX Z spin chain is shown 
to be given by power-functions of the distance between spin operators. This result is nothing 
but the theorem in section |l|. 

4 Concluding remarks 

We have analyzed the operator- valued Riemann-Hilbert problem associated with the generating 
functional of correlation functions and have obtained the large-m asymptotic form of a two- 
point correlation function (cr^crf) of the XX Z spin chain. This correlation functions is exactly 
expressed by power-functions. Our approach based on Riemann-Hilbert problems is useful for 
the analysis of long distance correlations. In a forthcoming publication we will determine the 
exponent of power, calculating the integrals ( |3.18 ) and ( p.19 ). The exponent is expected to 



depend on the anisotropy parameter rj and the external magnetic field h. 

As stated in section [I], a special correlation function FSFP is obtained from the generating 
functional by taking the limit a — > — do. However, in this limit, our asymptotic expression ( |3.23| ) 
loses its meaning. Similar situation arises from the asymptotic analysis of FSFP of the XXX 
spin chain (A = 1) 111]]. From the viewpoint of Riemann-Hilbert problem the evaluation of 
FSFP is more difficult than that of the generating functional. It is an open problem to compute 
the asymptotic form of FSFP for the XXX and XX Z spin chain. It is worth mentioning that 
in the free fermion point A = the corresponding FSFP has been shown to decay in Gaussian 
form |2). 

In this paper we have deformed an operator- valued Riemann-Hilbert problem for 4x4 matrix 
integral operators to scalar problems, which has been studied in detail |l(J. In general, it is hard 
to solve operator-valued Riemann-Hilbert problem: even in the case that the conjugation matrix 
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has no integral operator the solution can not be found in the closed form, and for operator-valued 
case its existence and uniqueness have not yet been proved. They are fundamental and unsolved 
questions. The investigation of operator-valued Riemann-Hilbert problem is important since it 
contains many applications to mathematics and physics. 

Acknowledgment 

The authors wish to thank Professor T. Deguchi and Dr T. Tsuchida for useful comments. 

References 

[1] V. E. Korepin, N. M. Bogoliubov, A. G. Izergin: "Quantum Inverse Scattering Method and 
Correlation Functions" , (Cambridge University Press, 1993, Cambridge). 

[2] F. H. L. Efiler, H. Frahm, A. G. Izergin and V. E. Korepin: Commun. Math. Phys. 174 
(1995) 191-214. 

[3] F. H. L. Efiler, H. Frahm, A. R. Its and V. E. Korepin: J. Phys. A 29 (1996) 5619-5626. 
[4] Y. Fujii and M. Wadati: J. Phys. Soc. Jpn 68 (1999) 3227-3235. 

[5] F. H. L. Efiler, H. Frahm, A. R. Its and V. E. Korepin: Nucl. Phys. B 446 (1995) 448-460. 
[6] Y. Fujii and M. Wadati: J. Phys. A 33 (2000) 1351-1361. 

[7] M. Jimbo M and T. Miwa: "Algebraic Analysis of Solvable Lattice Models", (CBMS Re- 
gional Conference Series in Mathematics 85, 1995, AMS). 

[8] N. Kitanine, J. M. Maillet and V. Terras: Nucl. Phys. B 567 (2000) 554-582. 

[9] F. H. L. Efiler, H. Frahm, A. R. Its and V. E. Korepin: J. Phys. A 30 (1997) 219-244. 

[10] N. I. Muskhelishvili: "Singular Integral Equations: Boundary Problems of Function Theory 
and Their Applications to Mathematical Physics", (Second edition, Dover Publications, 
1992, New York). 

[11] H. Frahm, A. R. Its and V. E. Korepin: "An operator-valued Riemann-Hilbert problem 
associated with the XXX model", (CRM Proc. Lecture Notes 9, AMS, 1996) 133-142. 

[12] P. Deift, A. R. Its and X. Zhou: Ann. Math. 146 (1997) 149-235. 



9 



